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Abstract This study investigated the systematics of the megophryid genus Scutiger from eastern and western Nepal 
using molecular and morphological data. Our results support two divergent lineages, one of which has nuptial spines on 
the dorsal surface of the first three fingers while the other has spines only on the dorsal surface of the first two fingers. 
The Ghunsa lineage from eastern Nepal shows significant morphological and molecular differences to other species of 
genus Scutiger and is here described as a new species. Based on the molecular analysis, the Muktinath lineage from 
western Nepal is confirmed to be Scutiger boulengeri and represents a species complex widespread throughout the 
Himalayan region. The newly described taxon is endemic to the eastern Himalayas and currently known only from the 
Ghunsa valley, Taplejung district, Nepal. 
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1. Introduction the world’s largest elevation gradient from 60 to 8848 

meters, within a 150 -200 kilometers (Dobremez 1976). 
The Eastern Himalayas is one of the biological hotspots Fast flowing rivers created sharp and deep river gorges in 
in Asia and a place of endemism where new species of this region promoting geographic isolation that hinders 
vertebrates and invertebrates have been discovered in ` gene flow, thus leading to speciation and endemism 
recent years (Agarwal et al., 2014; Khatiwada et al., (Agarwal et al., 2014; Che et al., 2010). Amphibians, 
2017; Khatiwada et al., 2015; Mahony et al., 2011). The ` being poor dispersers and generally philopatric (Smith 
Kangchenjunga Conservation Area (KCA) is a part ofthe ` And Green 2005), and may be expected to have increased 
Eastern Himalayas with high topographical, hydrological 
and climatic variability (WWF 2007) and is rich in 
floral and faunal diversity. Furthermore, this region is 


rates of genetic differentiation in this region compared to 
other vertebrate species. 
Members of the genus Scutiger Theobald, 1868 have 


situated at the transition zone of Palaearctic and Oriental been reported from southwestern China (including Tibet), 


zoogeographical realms (Olson et al., 2001) and therefore 
species in this region have a close affinity with both 
realms (Che et al., 2010). The Eastern Himalaya occupies 


northern Myanmar, northern India and Nepal (Fei et al., 
2012; Frost 2018). Twenty-one species of Scutiger are 
known, and among these are species in the mountains 
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(Table S1). Over the last 15 years several species in 
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described from the Eastern Himalayas and adjacent 
region (Jiang et al., 2016; Nguyen et al., 2013; Orlov 
et al., 2015; Sondhi and Ohler 2011; Yang et al., 2016), 
reflecting an increase in amphibian taxonomic diversity. 
Nevertheless, only one species of Scutiger has recently 
been described from this region (Jiang et al., 2016). In 
Nepal, two megophryid genera, Xenophrys Günther, 1864 
and Scutiger are known, with four species of Scutiger 
toads reported (Dubois 1974; 1987). The first record of a 
Nepalese Scutiger was Scutiger boulengeri, described as 
Cophophryne alticola from the Kharta Valley by Procter 
(Procter 1922) and also reported by others (Dubois 1974; 
1987; Nanhoe and Ouboter 1987) from Mukthinath. 
Scutiger sikimmensis (Blyth 1854) was first recorded 
from Langtang by Smith (1951), later by Dubois (1974) 
with new collections from the same locality. Dubois 
(1987) reported specimens from a series of Nepalese 
locations (Thankur, Ghorapani, Langtang, Namche Bazar 
and Jalijale). Nanhoe and Ouboter (1987) confirmed its 
presence in Ghorapani area of Annapurna Conservation 
Area. The specimens reported by Smith and Battersby 
(1953) from Maharigaon (western Nepal) were later 
allocated to Scutiger nyingchiensis, now S. occidentalis 
by Dubois (1987). Scutiger nepalensis was described 
by Dubois in 1974 based on specimens from Khaptad 
area, and has since been collected from Dharpatan 
Dhule in western Nepal (Dubois 1987);Nanhoe and 
Ouboter (1987) extended the distribution to Gurja Ghat. 
Recently, Hofmann et al., (2017) carried out a molecular 
phylogenetic analysis of genus Scutiger from all across its 
distributional range using both mitochondrial and nuclear 
markers. Their study provided updated genetic insight of 
genus Scutiger from Nepal and other areas, and further 
recognized Scutiger nepalensis and Scutiger sikkimensis 
from the Nepalese Himalaya as highly divergent. 

During field surveys in 2015 and 2016 in the 
Kanchenjunga Conservation Area (Taplejung district, 
eastern Nepal) and in the Annapurna Conservation Area 
(Mustang district, western Nepal), we collected a series 
of adult specimens and tadpoles that were assigned to 
the genus Scutiger based on the following characters: 
(1) enlarged paratoids; (2) granular dorsal surface; (3) 
absence of vomerine and maxillary dentition; (4) hidden 
tympanum; (5) finger and toes free of webbing; (6) 
absence of digital pads on fingers and toes; (7) fingers 
of breeding males with black nuptial spines and with 
pectoral and axillary glands. Detailed morphological and 
molecular analyses suggested that Ghunsa specimens are 
genetically and morphologically distinct from previously 
known species of genus the Scutiger. Thus, we provide a 


description of the Ghunsa population of Scutiger as a new 
species. 


2. Materials and Methods 


2.1. Sampling Fieldwork was carried out between May 
to September in 2015 and 2016 in eastern and central 
Nepal (Figure 1). Multiple survey methods were used; 
for example, visual encounter surveys, acoustic surveys, 
and leaf litter searches were made during the night 
(1900h—0000h) as described by Khatiwada et al., (2016) 
and Khatiwada et al., (2019). Specimens were collected 
by hand, euthanized using 20% Benzocaine gel, tissue 
samples were collected, and specimens later fixed in 4% 
formalin for 24 h and then preserved in 70% ethanol. 
Tadpoles and eggs at different developmental stages were 
also sampled from the same location as the type series. 
These samples were preserved in 70% alcohol. Tissue 
samples were taken by clipping the fourth toe from adults 
and small portion of tail muscle of tadpoles, which were 
then preserved to 95% ethanol for molecular analysis. All 
amphibians handling and processing were in accordance 
with the guidelines of the Department of National Park 
and Wildlife Conservation, Nepal and Animal Care and 
Use Committee of Chengdu Institute of Biology, CAS. 


2.2. Morphological analyses Morphological measure- 
ments were taken using digital calipers (to nearest 0.1 
mm). Measurement descriptions are provided in Table 
1. Preserved Scutiger specimens were also measured 
in the herpetological museum of Chengdu Institute of 
Biology (CIB), Chinese Academy of Sciences, Chengdu, 
China, and at the Collection des Reptiles et Amphibians 
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Figure 1 Study area map showing sampling sites in Nepal 
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Table 1 Morphological characters used and their measurement descriptions for adults and tadpoles. 


Abbreviation Morphology 


Measurement description 


Adult morphology 


SVL 


Snout vent length 


Distance from tip of snout to posterior edge of vent 


HL Head length Distance from angle of jaws and snout-tip 

HW Head width Measured at posterior angle of jaws 

SL Snout length From tip of snout to anterior corner of eye 

SN Nostril to the tip of the snout Distance from nostril to tip of snout 

ED Eye diameter Horizontal diameter of eye 

IND Internarial distance Minimum distance between external nares 

ENL Eyelid-naris distance Minimum distance between eyelid and rim of naris 

UEW Upper eyelid width Greatest width of upper eyelid 

IOW Inter orbital width Minimum distance between upper eyelids 

TYD Tympanum diameter Largest tympanum diameter 

TYE Tympanum orbit distance Distance from anterior border of tympanum to posterior orbital border 
FAL Forearm length Distance from elbow to base of outer metacarpal tubercle 

HAL Length of hand Measured from base of outer metacarpal tubercle to tip of third finger 
F1 First finger length Distance from base of second finger to tip of first finger 

F2 Second finger length Distance from base of first finger to tip of second finger 

F3 Third finger length Distance from base of second finger to tip of third finger 

F4 Fourth finger length Distance from base of third finger to tip of fourth finger 

THL Thigh length Distance from groin to knee 

SHL Shank length Distance from knee to heel 

TL Tarsus length Distance from heel to inner metatarsal tubercle 

FOL Foot length Distance from inner metatarsal tubercle to tip of fourth toe 

Tl First toe length Maximum length from base of first subarticular tubercle first toe tip 
T2 Second toe length Maximum length from base of first subarticular tubercle second toe tip 
T3 Third toe length Maximum length from base of first subarticular tubercle third toe tip 
T4 Fourth toe length Maximum length from base of first subarticular tubercle fourth toe tip 
T5 Fifth toe length Maximum length from base of first subarticular tubercle fifth toe tip 
Tadpole morphology 

TL Total length Distance from tip of snout to tip of tail 

BL Body length Distance from tip of snout to junction of posterior body and tail musculature 
BW Body width Highest body width 

BH Body height Highest body height 

HW Head width Measured at posterior angle of the jaws 

TAL Tail length Measured from junction of posterior body the tail musculature to tail tip 
TMH Tail muscle height at base of tail Maximum tail muscle height 

UTF Upper tail fin height Maximum height of upper fin, from upper margin of the tail musculature to upper margin of upper fin 
LTF Lower tail fin height Maximum height of the lower fin, from lower margin of lower fin to lower margin of tail musculature 
TMW Tail muscle width at base of tail Maximum tail muscle width 

MTH Maximum height of tail Maximum height of the tail 

ENL Eyelid-naris distance Minimum distance between eyelid and rim of naris 

SN Nostril to the tip of the snout Distance from the nostril to the tip of the snout 

SSD Snout-spiracle distance Distance from tip of snout to posterior margin of spiracle 

IND Internarial distance Minimum distance between external nares 

IOW Inter orbital width Minimum distance between upper eyelids 

ED Eye diameter Horizontal diameter of eye 


of Muséum national d’Histoire naturelle (MNHN), 
Paris, France. A total of 83 individuals (51 males and 32 


females) of Scutiger from nine species were examined 
(Table S2). 


2.3. Statistical analysis A one-way Analysis of Variance 
ANOVA was used to test for differences among the 


new species, A. boulengeri, and S. nepalensis, followed 
by Tukey’s HSD post hoc test for pairwise differences. 
Prior to Principal Component Analysis (PCA), all 
morphometric data were regressed against SVL and the 
regression residual of each morphological trait was used 
as input data for the PCA. Multivariate Analysis of 
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Variance (MANOVA) was used to test the morphometric 
differences between species by using scores from the 
principal components (those with eigen-value > 1) as 
dependent variables and species as independent variables 
(groups). Statistical analysis was performed in R package 
vegan (Oksanen et al., 2016) in R-software version 3.3.0 
(R Development Core Team 2016). 


2.4. Molecular systematics analyses The phylogenetic 
position of the Ghunsa population of Scutiger was 
determined in relation to other Scutiger species through 
analysis of DNA sequence data. Total genomic extraction 
was carried out on clipped toe or muscle preserved in 
95% ethanol using the DNeasy Tissue Kit (QIAGEN). 
The phylogenetic position of the Ghunsa population of 
Scutiger was determined in relation to other Scutiger 
species through analysis of DNA sequence data. The 
mitochondrial gene 16S ribosomal RNA gene (16S 
rRNA), complete cytochrome C oxidase | gene (COD and 
cytochrome b (cytb) from all samples were sequenced and 
amplified using primers as detailed by Che er al., (2012). 
The polymerase chain reaction (PCR) and thermocycling 
conditions were set as described by Che et al., (2012). The 
amplified PCR products were purified using Qiagen PCR 
purification kit and sequenced were obtained from an ABI 
3100 automated sequencer (Applied Biosystems). Newly 
obtained sequences were deposited in GenBank under the 
accession numbers (16s: MK950885—-MK950907, COI: 
MK970591—MK970613, Cytb: MK970614—MK97063 1). 
Available nucleotide sequences of Scutiger species were 
downloaded from the NCBI Genbank database and a list 
of species included in our analysis is provided in Table 
2. ClustalW built into BIOEDIT Version 7.1.9 was used 
to align the nucleotide sequences (Thompson et al., 
1997) using the default parameters. Alignments were also 
checked and manually edited as needed. 

Maximum likelihood (ML) and Bayesian inference 
(BI) analyses were employed to reconstruct the 
phylogenetic relationships among the taxa based on 
haplotypes of the concatenated sequences of three genes 
(16S, COI and cytb), and COI, and 16S individually. 
ML analysis was conducted with the rapid bootstrapping 
algorithm using the program RAxML v8.00 (Stamatakis 
2014) on the CIPRES Science Gateway server v3.2 
(Miller et al., 2010). Nodal support for ML was assessed 
with 1000 rapid bootstrap replicates (BS). 

PartitionFinder v1.1.0 (Lanfear er al., 2017) was used 
to find the optimal evolutionary models for Bayesian 
inference analysis based on Bayesian Information 
Criterion (BIC) values (Darriba et al.,, 2012). For 
the concatenated mtDNA, seven potential data block 


partitions were defined, one block for 16s rRNA gene and 
six blocks representing the three codon positions of each 
of the two protein-coding genes (CYTB and COI). BI 
analyses were conducted using MrBayes 3.1.2 (Ronquist 
and Huelsenbeck 2003). Two independent runs were 
initiated, each with four simultaneous Markov Chain 
Monte Carlo (MCMC) chains, for 20 million generations 
and sampled every 1000 generations. The adequacy 
of convergence of chains and burn-in period was 
determined by examining plots of log-likelihood scores 
and the standard deviation of split frequencies. , The 
first 25% generations were discarded as burn-in, and the 
remaining samples were used to create a 50% majority- 
rule consensus tree and estimate Bayesian posterior 
probabilities (BPP). The graphical viewer Figtree 
(Rambaut 2007) was used to edit the output of RAxML 
and MrBayes analyses. Uncorrected k2p pair-wise 
distances within and across genus for COI gene sequences 
were calculated with MEGA7 using the complete deletion 
option (Kumar et al., 2016). 


3. Results 


Phylogenetic relationships of mtDNA haplotypes The 
aligned dataset of concatenated DNA sequences of the 
mitochondrial genes (16S, COI and Cyt-b) had a total 
length of 1939 bp. The maximum likelihood (ML) and 
Bayesian inference (BI) phylogenetic trees constructed 
from concatenated DNA sequences produced similar 
topologies with strong node support (Figure 2). The 
phylogenetic analysis indicates that the population of 
Scutiger from Ghunsa, (Kanchenjunga Conservation 
Area) is weakly supported as the sister taxon to a clade 
composed of S. nepalensis and S. sikkimensis. Several S. 
boulengeri were identified and it was misidentified as A. 
mammutus by Hoffman et al., (2017). 

Based on the COI gene only, the Scutiger population 
from Muktinath (Mustang district, Annapurna 
Conservation Area, western Nepal) was genetically 
similar to Scutiger boulengeri from nearby neotype 
locality (Zhiduo and Zhadao counties, Qinghai, China) 
and clustered together in the phylogenetic tree (Figure 
S1). Furthermore, the uncorrected COI genetic distances 
among the Ghunsa population of Scutiger from the 
Kanchenjunga Conservation Area were <0.01%. The 
distances to the closet congener, S. nepalensis were 
7.4%—-11.8%. The genetic divergences between the 
Ghunsa population and other two known species of 
Scutiger in Nepal, S. sikkimmensis from Arun river basin, 
eastern Nepal, see in Hofmann er al., (2017) and S. cf. 
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(JRK193) Ghunsa. 
(JRK194) Ghunsa, 
(JRK195) Ghunsa, 
(JRK197) Ghunsa, 
(JRK199) Ghunsa, 
(JRK206) Ghunsa, 
(JRK202) Ghunsa, 
(JRK203) Ghunsa. 
(JRK196) Ghunsa. 
(JRK198) Ghunsa, 
(JRK200) Ghunsa. 
(JRK201) Ghunsa, 
(JRK204) Ghunsa, 
(JRK205) Ghunsa, 
(JRK207) Ghunsa, 


Scutiger ghunsa sp. nov. Taplejung, Nepal 


Scutiger ghunsa sp. nov. Taplejung, Nepal 


Scutiger ghunsa sp. nov. Taplejung, Nepal 


Scutiger ghunsa sp. nov. Taplejung, Nepal 


Scutiger ghunsa sp. nov. Taplejung, Nepal 


Scutiger ghunsa sp. nov. Taplejung, Nepal 


Scutiger ghunsa sp. nov. Taplejung, Nepal 
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Scutiger ghunsa sp. nov. Taplejung. Nepal 
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Scutiger ghunsa sp. nov. Taplejung, Nepal 


Scutiger ghunsa sp. nov. epal 


Scutiger ghunsa sp. nov. Taplejung, Nepal 


Scutiger ghunsa sp. nov. epal 


Scutiger ghunsa sp. nov. Taplejung, Nepal 


Scutiger ghunsa sp. nov. 
S. nepalensis, Kanjiroba, Jumla, Nepal 

S. nepalensis, Saipal-Juphal, Nepal 

S. nepalensis, Dhauli lake, Jumla, Nepal 

S. nepalensis, Dhaulung-North, Jumla, Nepal 

S. nepalensis, Dhaulung-South, Jumla, Nepal 

S. nepalensis, Saipal-chainpur, Jumla, Nepal 

S. sikkimensis, Nepal, Arun, Kongma Danda, Nepal 

S. sikkimensis, Taksindo, Dudhkoshi river-East, Nepal 
S. sikkimensis, Surkie La, Dudhkoshi river-West, Nepal 

S. chintingensis, Tianquan, Sichuan, China 

S. occidentalis, Deosai Plains, Pakistan 


0.64/16 


0.61/19 


S. boulengeri, Damxung, Tibet, China 
S. boulengeri, Lasha, Tibet, China 

S. boulengeri, Zhiduo, Qinghai, China 
S. boulengeri, Tageija, Tibet, China 

S. boulengeri, Garze, Sichuan, China 

S. boulengeri, Nyingchi, Tibet, China 


S. boulengeri, Gansu, China 
Orelalax major 


Figure 2 Bayesian inference (BI) tree based on concatenated DNA sequences of 16S, COI and Cytb. The values on branches are Bayesian 
posterior probabilities (bpp) followed by ML bootstrap supports (bs). 


boulengeri from Mustang, Annapurna Conservation Area, 
western, Nepal were 10.10% to 11.4% and 10.10% to 
10.3% respectively. 

Within S. nepalensis (the clade of six specimens 
from Jumla district-Hofmann et al., [2017]) and the 
sample from Saipal mountain, Chainpur, Jumla district, 
Nepal showed 5.5-5.8% genetic divergence from other 
individuals of this group. Within S. sikkimensis group, the 
sample from Yadong, Tibet, China (see details in Jiang 


et al., [2016]) showed 9.9% divergence from populations 
from Kongma Danda, Arun River, Taksindo, Dudhkoshi 
River (east), and Surkie La, (see in Hofmann et al., 
[2017]). Furthermore, within samples of S. sikkimensis 
in Nepal, Surkie La, Dudhkoshi river (west) population 
showed 7.1% genetic divergence from the Kongma 
Dandaand Taksindo, populations (Figure S1). 

Based on 16S sequences alone, S. sikkimensis from 
Sikkim, India showed lower genetic distances (0.00- 
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Table 2 Samples used in molecular analysis with GenBank accession numbers. 


GenBank GenBank GenBank 
Species Voucher number Locality accession accession accession Source 


No. (COI) No. (cytb) number (16s) 
Ingroup 
S. ghunsa sp. nov. JRK2015-193 Ghunsa, Taplejung district, Nepal MK970591 MK970614 MK950885 This study 
S. ghunsa sp. nov. JRK2015-194 Ghunsa, Taplejung district, Nepal MK970592 MK970615 MK950886_ This study 
S. ghunsa sp. nov. JRK2015-195 Ghunsa, Taplejung district, Nepal MK970593 MK970616 MK950887_ This study 
S. ghunsa sp. nov. JRK2015-196 Ghunsa, Taplejung district, Nepal MK970594 MK970617 MK950888 This study 
S. ghunsa sp. nov. JRK2015-197 Ghunsa, Taplejung district, Nepal MK970595 MK970618 MK950889 This study 
S. ghunsa sp. nov. JRK2015-198 Ghunsa, Taplejung district, Nepal MK970596 MK970619 MK950890 This study 
S. ghunsa sp. nov. JRK2015-199 Ghunsa, Taplejung district, Nepal MK970597 MK970620 MK950891 This study 
S. ghunsa sp. nov. JRK2015-200 Ghunsa, Taplejung district, Nepal MK970598 MK970621 MK950892_ This study 
S. ghunsa sp. nov. JRK2015-201 Ghunsa, Taplejung district, Nepal MK970599 MK970622 MK950893 This study 
S. ghunsa sp. nov. JRK2015-202 Ghunsa, Taplejung district, Nepal MK970600 MK970623 MK950894_ This study 
S. ghunsa sp. nov. JRK2015-203 Ghunsa, Taplejung district, Nepal MK970601 MK970624 MK950895_ This study 
S. ghunsa sp. nov. JRK2015-204 Ghunsa, Taplejung district, Nepal MK970602 MK970625 MK950896_ This study 
S. ghunsa sp. nov. JRK2015-205 Ghunsa, Taplejung district, Nepal MK970603 MK970626 MK950897_ This study 
S. ghunsa sp. nov. JRK2015-206 Ghunsa, Taplejung district, Nepal MK970604 MK970627 MK950898 This study 
S. ghunsa sp. nov. JRK2015-207 Ghunsa, Taplejung district, Nepal MK970605 MK970628 MK950899 This study 
S. boulengeri JRK2016-211 Muktinath, Mustang district, Nepal MK970606 - MK950900 This study 
S. boulengeri JRK2016-212 Muktinath, Mustang district, Nepal MK970607 - MK950901 This study 
S. boulengeri JRK2016-213 Muktinath, Mustang district, Nepal MK970608 - MK950902 This study 
S. boulengeri JRK2016-214 Muktinath, Mustang district, Nepal MK970609 - MK950903 This study 
S. boulengeri JRK2016-215 Muktinath, Mustang district, Nepal MK970610 - MK950904 This study 
S. boulengeri JRK2018-01 Zhiduo, Qinghai, China MK970611 MK970629 MK950905_ This study 
S. boulengeri JRK2018-02 Zhiduo, Qinghai, China MK970612 MK970630 MK950906_ This study 
S. boulengeri JRK2018-03 Zhaduo, Qinghai, China MK970613 MK970631 MK950907_ This study 
S. cf. boulengeri AI-AL a a, Tibet Autonomous Region, 4319879 KY310922 KY310760 Hofmann et al. (2017) 
S. cf. boulengeri A3-AL Damne, Tibet Autonomous KY310872 KY310924 KY310762 Hofmann et al. (2017) 
Region, China 
S. cf. boulengeri JS1507 Cl Se Tibet Autonomous Region, KY310875 KY310927 KY310765 Hofmann ef al. (2017) 
S. cf. boulengeri KQ2 2014 Garze Sichuan, China KY310860 KY310912 KY310750 Hofmann er al. (2017) 
S. cf. boulengeri GAs snipes er Tibet Autonomous KY310877 KY310929 KY310768 Hofmann etal. (2017) 
Region, China 
S. cf. boulengeri jonel Jone County, Gansu, China KJ082073 KJ082065 EF39727 Hofmann et al. (2017) 
S. cf. boulengeri B07 Chopta Valley, Sikkim, India - - KR781481 Subba et al. (2015) 
S. nepalensis NME_A2014/13 Kanjiroba, Jumla, Nepal KY310936 KY310885 KY310776 Hofmann et al. (2017) 
A. nepalensis NME_A2018/13 Chainpur, Jumla, Nepal KY310937 KY310886 KY310777 Hofmann et al. (2017) 
S. nepalensis A8-2012 Dhaulgiri north, Dolpa, Nepal KY310938 KY310887 KY310778 Hofmann et al. (2017) 
S. nepalensis A14-2012 Dhaulagiri, south, Dolpa, Nepal KY310944 KY310893 KY310784 Hofmann et al. (2017) 
S. nepalensis MS_Scut2 Juphal, Dolpa, Nepal KY310947 KY310896 KY310787 Hofmann et al. (2017) 
S. nepalensis MS _Scut5 Dhauli Lake Nepal KY310948 KY310897 KY310788 Hofmann et al. (2017) 
S. sikkimensis JS140524 SE Danda, Sankhuwasava, ` Kaz10e? KY310959 KY310798 Hofmann et al. (2017) 
S. sikkimensis JS140528 Surkie La, Solukhumbu, Nepal KY310904 KY310963 KY310801 Hofmann et al. (2017) 
S. sikkimensis Ne11-2013 Taksindo, Solukhumbu, Nepal KY310909 KY310969 KY310806 Hofmann etal. (2017) 
S. sikkimensis KIZ011127 aia Tibet Autonomous erguenen 2: R Jiang et al. (2016) 
S. sikkimensis B07 Chopta valley, Sikkim, India - - KR781486 Subba et al. (2015) 
S. occidentalis MS _PK2 Deosai Plains, Pakistan KY310899 KY310954 KY310794 Hofmann et al. (2017) 
S. chintingensis LC141 Tianquan, Sichuan, China KY310878 KY310930 KY310769 Hofmann et al. (2017) 
S. glandulatus SC1_2014 Garze, Sichuan, China KY310879 KY310934 KY310770 Hofmann et al. (2017) 
S. liupanensis FWRI11 Liupan Mountains, China KC140483 JX533792 KX352261 Hofmann et al. (2017) 
S. ningshanensis bys1 Taibai, Shaanxi, China KF75739 KF757366 KX619450 Hofmann et al. (2017) 
A. tuberculatus CIBXM988 Sichuan, China - FJ945493 EF397278 Hofmann et al. (2017) 


S. gongshanensis CIB20070717001 Gongshan, Yunnan, China KU243062 - EF397278 Jiang et al. (2016) 
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(Continued Table 2) 
GenBank GenBank GenBank 
Species Voucher number Locality accession accession accession Source 
No. (COI) No. (cytb) No. (16s) 
S. jiulongensis KIZ045055 Ganzi, Sichuan, China KU243066 - - Jiang et al. (2016) 
n ee Medog, Nyingchi, Tibet . 

S. nyingchiensis KIZ017459 Autonomous Region, China KU243056 - - Jiang et al. (2016) 

or Medog, Nyingchi, Tibet : 
S. spinosus KIZ011114 Autonomous Region, China KU243053 - - Jiang et al. (2016) 

f , Medog, Nyingchi, Tibet ? 
S. wuguanfui KIZ011101 Autonomous Region, China KU243060 - - Jiang et al. (2016) 
S. muliensis IOZCAS3638 Yanyuan, Sichuan, China - - EF397277 Fu et al. (2007) 
Outgroup 
Oreolalax major CIB-LX548 Luding, Sichuan, China KU127230 KU127230 - Liang et al. (2016) 
Xenophrys omeimontis KIZ025765 Emei Shan, Sichuan, China KX812136 KR087950 - Fu et al. (2007) 
Oreolalax xiangchengensis ` - Li Jiang, Yunnan, China - - NC032349 - 
Xenophrys jingdongensis 1999.5688 Sapa, Lao Cai, Vietnam KR087947 - - Grosjean et al. (2015) 
Megophrys maosonensis KIZ016045 Xiaoqiaogo Nature Reserve, KX812080 - - Chen et al. (2017) 

Yunnan, China 

Duttaphrynus melanostictus  - - AY458592 AY458592 AY458592 Zhang etal. (2005) 


0.20%) from of the Kongma Danda and Taksindo) 
populations (Figure S2). Scutiger boulengeri from 
Muktinath, Mustang, Nepal were genetically very 
similar with the specimens from Zhiduo and Zhadao, 
Qinghai, China and Tageija, Damxung and Lasha, Tibet, 
China (Figure S2) with genetic distances of only 0.1 to 
0.3%. Scutiger boulengeri from Sikkim, India were also 
genetically very similar with Scutiger from Muktinath, 
Nepal; the genetic distance between them was 0.4% 


Morphometric analysis Morphometric comparisons 
were used to verify the distinctiveness of the lineages 
identified by molecular studies. Based on phylogenetic 
analysis, the Scutiger population from Ghunsa was 
closely related to S. nepalensis. Therefore, quantitative 
and multivariate analysis was performed from three 
species of Scutiger (Scutiger Ghunsa population, S. 
nepalensis and S. boulengeri) known to occur in Nepal. 
In Scutiger Ghunsa population, females were significantly 
larger than males in SVL (t = —4.50, p = 0.006, df = 8). 
Due to lack of females of S. nepalensis, we only used 
males for the further morphological comparisons. One- 
way ANOVA showed that mean SVL varied significantly 
among the three species, ranging from 33.8 to 76.0 mm 
among males (one-way ANOVA, p < 0.001; F = 177.3). 
Further, post hoc analyses revealed that the mean SVL 
was significantly different among the three species 
(Figure 3). 

The Principal Component Analysis extracted four 
components with eigenvalues > 1.0 which accounted 
for 91.7% of the variation in the dataset. The first three 
components explained 85.6% of total variation. In PCA 
plots, the first two principal components (PC1 vs PC2; 


Figure 4A) and the first vs. the third principal components 
(PC1 vs PC3; Figure 4B) for males separated the three 
species from each other. Species variation along PC 
describes essentially the variation in morphological traits 
with longer length of the nostril to eye distance, inter 
orbital width, inter narial width, length of hand, 2nd 
finger, 4th finger, length of tarsus, length of meta-tarsus, 
Ist toe, 2nd toe, 3rd toe, 4th toe, 5th toe having positive 
scores whereas, species with shorter morphological 
traits have negative scores along the PC2 (Table 3). All 
the three Scutiger species were significantly different in 
morphometric shape (MANOVA, Wilks’ à = 0.035, F= 
8.57, p < 0.001). The difference was in accordance with 
the PCA results. 


Systematics These results confirmed that there is 
substantial genetic divergence between the specimens 
from Ghunsa and other known Scutiger species, and 
support that this population is a distinctly evolving 
lineage, representing an undescribed species. This 
population also presents distinct morphological 
characteristics that are not observed in the closely related 
species S. nepalensis, which has smaller nostril-to-eye 
distance, inter-orbital width, inter-narial width, shorter 
length of hand, 2nd finger, 4th finger, length of tarsus, 
length of meta-tarsus, 1st toe, 2nd toe, 3rd toe, 4th toe and 
5th toe. Based on molecular and morphological evidence, 
we hereby describe the specimens collected from Ghunsa, 
Taplejung district, Kangchenjungha Conservation Area, 
Nepal as a new species. 


Species Account: Scutiger ghunsa sp. nov. (Figure 5) 
Local Name: Ghunsa alpine toad 
Nepali Name: Ghunsa Khasre Bhyaguto 
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Figure 3 Body size (mean+SD) variation among male Scutiger 
species. The bar represents the mean SVL and error denotes the 
standard deviation. The statistics derived from a one-way ANOVA 
with post hoc Tukey test. The different letters above the bar signify 
the significant difference at 0.05. 


Holotype: NHM-TU-17A-0116 (Figure 5), an adult male 
collected from Ghunsa Khola (Khola means “stream” in 
Nepali) Ghunsa, Taplejung district, Nepal, 27.66185°N, 
87.93168 °E, elevation 3475 m a.s.l. collected by Janak 
Raj Khatiwada between 19:30 and 20:00h on 28" May 
2015 and deposited in the collection of the Natural 
History Museum, Tribhuvan University, Soyambhu, 
Kathmandu, Nepal. 


Suggested common name: Ghunsa High Altitude Toad 


Etymology: The species name is derived from the name 
of the type locality. It is hoped that this name will raise 
awareness to the local people about this new species and 
will further help to conserve the flora and fauna in this 
area. The specific epithet is a noun in apposition and is 
invariable in gender. 


Paratypes: NHM-TU-17A-0117 (adult male) and NHM- 
TU-17A-0118 (adult female, Figure 6) from the same 
location of as the holotype. The specimens are deposited 
in the collection of the Natural History Museum, 
Tribhuvan University, Soyambhu, Kathmandu, Nepal. 


Diagnosis: Scutiger ghunsa sp. nov. is assigned to the 
genus Scutiger, based on both molecular phylogenetic 
data and the following morphological characters: (1) 
enlarged paratoids; (2) dorsal surface granular in adult 
males; (3) vomerine and maxillary dentation absent; (4) 
hidden tympanum; (5) finger and toes free of webbing; (6) 
absence of digital disk on fingers and toes; (7) inner three 
fingers of breeding males with black nuptial spines in 
adult males; (8) breeding males with pectoral and axillary 
glands. 

Scutiger ghunsa sp. nov is distinguished from all 


other known species of Scutiger by having the following 
morphological characters: (1) adult with medium body 
size, SVL 42.1-47.2 mm in males and 50.2—53.9 mm in 
females; (2) head wider than longer; (3) snout shorter than 
20-40% of head length; (4) eye shorter than 20% of head 
length; (5) interorbital distance longer than internarial 
distance; (6) maxillary teeth absent; (7) vocal sac absent; 
(8) odontoids present on the lower mandibles; (9) hidden 
tympanum; (10) length of arm 20—40% of SVL; (11) hand 
length shorter than 20% of SVL (12) breeding males with 
pectoral and axillary glands; (13) inner pectoral gland 
almost twice the size of outer gland in breeding males; 
(14) nuptial spines on dorsal surface of first, second 
and lateral surface on third fingers in breeding males; 
(15) finger and toes free of webbing; (16) subarticular 
tubercles absent (17) dorsal surface of body and limbs 
granular; (18) belly of breeding males without spines. For 
detailed morphological comparisons, see Table 4. 


Description of holotype (all measurements in mm) 
Moderate body size (SVL 43.6); head wider than long 
(HW 16.4, HL 11.4, HW:HL 143 %); snout short and 


Table 3 Summery table of PCA statistics showing: character 
loadings, Eigenvalues, and percentage of explained variance for 
Principal Components (PC) 1—4. 


EE PC 1 PC2 PC3  PC4 
HL —0.069 0.868 -0.382 0.103 
HW -0.464 0.649 0.506 0.102 
SL -0.322 0.889 0.188 mum 
SN 0.753 -0.296 0.175 0.527 
ED 0.971 0.092 -0.063 —0.147 
IND 0.782 ` mune 0.575 0.162 
ENL 0.991 -0.059 0.050 0.005 
UEW 0.261 0.218 0.209 0.894 
IOW 0.399 0.529 0.371 —0.168 
TYD 0.993 0.049 0.024 0.027 
TYE 0.448 0.645 -0.151 017 
FAL 0.979 0.121 0.011 -0.036 
HAL 0.577 0.445 0.427  —0.098 
FI 0.960 0.128 0.090 -0.173 
F2 -0.329 0.409 -0.77 0.337 
F3 0.778 0.46 0.0 0.15 

F4 0.942 0.064 0.075 0.066 
THL 0.857 0.307 -0.292 -0.05 

SHL 0.993 -0.014 0.058 0.035 
TL 0.993 0.008 -0.073 0.012 
FOL 0.976 -0.02 -0.143 0.02 

Tl 0.877 0.033 0.434 0.09 

T2 0.984 -0.031 0.078 -0.023 
ne vee 61.10 76.41 85.59 91.72 
Eigenvalue 14.01 3.52 2.11 1.41 
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Figure 4 Principal Component Analysis (PCA) on male specimens of Scutiger from Nepal based on 24 morphometric traits. Prior to 
Principal Component Analysis (PCA), all morphometric data were regressed against SVL and regression residual of each morphological trait 


was used in a PCA. 


round (SL:HL 41.6%); large paratoid glands on either side 
of head between posterior corner of eye to armpit (axills) 
(7.1 mm); canthus rostralis distinct; nostril dorsolateral, 
just below canthal ridge; loreal region slightly concave, 
eyes large and convex (ED 6.3); eye diameter almost 
similar to snout length (ED:SL 99 %); tympanum and 
supratympanic fold absent; interorbital space flat, 
interorbital distance (IOD 5.3), greater than width of the 
upper eyelid (UEW 4.3); internarial surface flat (ENL 
4.0); choanae small; vomerine teeth and ridges absent; 
maxillary teeth absent; tongue oval, deeply notched 
posteriorly, lingual papillae absent; openings of vocal sacs 
absent. 

Arm robust, short (FAL:SVL 23%) and shorter 
than hand (HAL:SVL 54%); relative finger lengths 
I<II<IV<IU; third finger (F3 5.9) shorter than arm (FAL 
10); fingers tip round, not dilated; fingers without distinct 
lateral fringes, webbing absent; subarticular tubercles 
indistinct; flat, large and ovoid metacarpal tubercles, inner 
tubercle almost equal in size to outer tubercle; nuptial 
spines on dorsal and lateral surface of first and second 
fingers, but only on inner side of third finger. 

Hindlimbs powerful and long, tibiotarsal articulation 
reaches tympanum when hindlimb pressed parallel to 
body; shank (SHL 17.5) almost equal to thigh (THL 
17.0) but tarsus short (TL 11.1), foot length longer (FOL 
19.2) than shank and thigh; toes thin and long, relative 
length I<II<III<V<IV (fourth toe of right foot cut for 
molecular analysis); toe tips round; weakly webbed; 
subarticular tubercles indistinct; inner metatarsal tubercle 
flat and oval; outer metatarsal tubercle poorly developed; 


supernumerary and plantar tubercles absent. 


Skin texture: Skin of dorsum rough with scattered warts 
bearing black spines, density of warts and granules 
increasing towards vent; dorsal surface of snout and 
interorbital region smooth; few scattered tubercles on 
upper and lower mandibles; dorsal surface of forelimbs 
and hindlimbs with smaller tubercles than back; flanks 
with larger white warts or granules; side of shank with 
white warts or glandular tubercles without spines; a pair 
of pectoral and axillary glands present on chest, pectoral 
glands two times longer than axillary glands, glands 
covered by minute, dense black spines; vent surrounded 
by granules. 


Colouration in life: An olive brown triangular spot on 
snout and interorbital region; posterior surface of head 
and dorsal body light brown; alternating olive brown and 
dark brown bands on upper lip; iris golden; forelimbs 
and hind limbs olive brown; irregular black spots on both 
limbs including fingers and toes; flanks light brown and 
gradually fading into creamy yellow ventrally. Ventral 
surface of head light brown; chest and abdomen creamy 
white with irregular light gray-brown lines; ventral 
surface of limbs brown with small irregular creamy white 
spots. 


Colouration in preservation: In ethanol, olive brown 
triangular spot on snout and interorbital region fading 
to white; dorsum body, head and upper limbs brownish- 
black; ventral surface of belly creamy white; limbs light 
brownish-black; black spines more distinct on shank and 
thigh; and iris black. 
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Figure 5 Holotype of Scutiger ghunsa sp. nov. (NHM-TU-17A-0116, adult male). in life: A. dorsolateral view, B. ventral view. In 
preservative: C. dorsal view, D. ventral view, E. lateral view, F. ventral view of chest, G. Dorsal surface of hand, H. ventral surface of foot. 


Janak Raj KHATIWADA et al. 
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Figure 6 Paratype of Scutiger ghunsa sp. nov. NHM-TU-17A-0117, adult female, SVL 47.8 mm. A. dorsal view; B. ventral view; C. ventral 


surface of hand; and D. ventral surface of foot. 


Variation and sexual dimorphism: Morphometric 
variation of type specimens is presented in Table 5. The 
female were significantly larger in SVL (t = -4.50, p = 
0.006, df= 6), head width (t = 2.63, p = 0.039, df= 6), 
inter narial width (t = 2.74, p = 0.034, df= 6), 3rd finger 
(t = 3.18, p = 0.019, df= 6), and 4th finger (t = 3.21, p= 
0.016, df = 6) than males. Males with nuptial spines on 
the dorsal surface of first three fingers (absent in females); 
a pair of pectoral and axillary glands present on chest 
(absent in females). Females with smooth skin, body 
surface devoid of spines; warts flatter on dorsal body 
surface than in males. 


Eggs 
A total of 14 eggs were collected while being guarded by 


a male under a rock. The eggs were rounded in shape and 
ranged from 3.4 to 4.2mm (Mean 3.9 + SD 0.28, N = 14). 
The vegetal and animal poles were creamy in colour. 


Larvae (all measurements in mm) 

Thirteen tadpoles (identity confirmed by molecular 
analysis) from different stages were examined (Gosner 
1960). The following description is based on Gosner stage 
27 tadpole of Scutiger ghunsa sp. nov. In dorsal view, 
body oval, posterior part wider than the anterior part, total 
body length (N = 5, TL = 27.6—33.1), body length 33% of 
TL and tail length 66% of TL. Head oval and broad, snout 
rounded, eyes dorsolaterally set, nostrils oval, closer to 
snout (ENL 18% of BL, SN 11% of BL, SSD 56% of BL, 
IND 22% of BL and IOW 29% of BL) than to eye. Oral 
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Table 4 Mensural and meristic data for adults and tadpole of Scutiger ghunsa sp. nov. compared with congeneric species. Morphological 
characters: 1. Male SVL (ranges with total number examined); 2. Female SVL (ranges with total number examined); 3. Maxillary teeth: 
(1) present, (0) absent; 4. Relative length of snout (SL/HL*100): (1) shorter than 20%, (2) between 20%-40%, (3) longer than 40% of HL; 
5. Relative eye size (ED/HL*100): (1) shorter than 20%, (2) between 20%—40%, (3) longer than 40% of HL; 6. Interorbital distance (IOW/ 
IND*100: (1) longer, (2) equal, (3) shorter than IND; 7. Length of arm (LA/SVL*100): (1) shorter than 20%, (2) between 20%-40%, (3) 
longer than 40% of SVL; 8. Length of hand (LA/SVL*100): (1) shorter than 20%, (2) between 20%-40%, (3) longer than 40% of SVL; 9. 
Number of fingers with nuptial spines: (1) presence in first three fingers (I, II, III), (2) first two fingers (I, II); 10. Inner metatarsal tubercle: (1) 
longer than first toe, (2) equal, (3) 50% less or shorter; 11. Subarticular tubercles: (1) present, (0) absent; 12. Toes webbing: (1) webbed, (2) 
weakly webbed, (0) absent; 13. Number of pectoral and axillary glands: (1) one, (2) two; 14: Size of pectoral and axillary glands: (1) pectoral 
and axillary glands same size, (2) pectoral glands slightly larger than axillary glands, (3) pectoral glands two time or even larger than axillary 
glands; 15. Spines on belly in males: (1) present, (0) absent; 16. Vocal sacs: (1) present (0) absent; 17. Keratodont row formulae in tadpole. 


Scutiger 1 2 3 4 5 6 7 8 9 10 1 12 13 14 15 16 17 Sources 
species 
ghunsa so se iad 0 2 3 1 2 3: L003 0 0 2 3 0 O 1:2+2/2+2:1This study 
(n=5) (n= 3) 
adungensis ZEN = 1 2 2 - 0 1 - Dubois (1979) 
(n= 2) 
n 64.9 This study; 
bhutanensis @=1) — 0 2 2 3 b3 LW 3 © 2 2 1 0 8 - Delorme and Dubois (2001) 
48.7 42.3 1 : 5+5 (6+6) / This study; 


© 
j 
N 
= 
= 
w 


LEM 3 0 1 2 1 1 


© 


boulengeri 


(n=5) (n=3) 5+5 (6+6): 1 Fei et al. (2012) 


; 68.0-80.0 68.0-80.0 : 
brevipes (n= 10) (n=9) 0 - - - - LD - - 2 2 3 0 0 - Fei et al. (2012) 


Say a, 38.1-44.8 47.5-50.3 aus _, This study; 
chintingensis (n=5) (n=4) 1 2 3 1 1 3 LIM 3 0 1 2 3 0 O 1:34+3/242: l Fei et al. (2012) 
55.7-80.8 62.8-82.5 wads _ , This study; 
glandulatus (a=5) Grady 0 2 3 1 2 3 II 3 1 2 2 3 0 O 1:444/414; l Fei et al. (2012) 


47.0-57.0 049.0—60.5 


gongshanensis =) (n= 2) 1 - - LD - 0 0 1 — O O 1:3+3/(+3) :1Fei et al. (2012) 
ai 67.4-81.5 1:3+3/4+4 |... 
jiulongensis (n= 20) — 0 - - - LU - 1 2 2 300 (343) :1 Fei et al. (2012) 


; ; 40.6-48.0 52.0-59.5 . . : 
liupanensis (n=20) (n=3) 1 - - - - -LIMH - 0 2 2 1 1 O 1:6+6/6+6: 1Fei et al. (2012) 


65.4 


maculatus get 69.0(n=1) 1 LIBH - 0 1 2 2 0 0 - Fei et al. (2012) 
57.5-65.0 67.3-75.2 1:4+4/5+5 This study; F 
mammatus (n=3) (n=5) 0 2 3 1 2 3 Ll 3 1 1 1 - 0 0 (4+4): 1 ei et al. (2012) 
aS 68.2-80.0 60.1—-67.5 ` _ , This study; 
muliensis (n= 10) (n= 10) 0 - - - - — LH - - 2 1 - O0 Q 1:414/4+4: l Fei et al. (2012) 
: 69.7-76.0 60.4-66.8 F 
nepalensis (n=4) (n=2) 7 2 2 3 2 2,0,M - - - 2 1 0 0 - This study 
i . 51.0 41.0 . P : 
ningshanensis =) =n 1 - - - - -LILII - - 2 2 1 1 O 1:5+5/5+5: 1Fei et al. (2012) 
, .  , 40.8-55.8 54.9-74.2 ` _ , This study; 
nyingchiensis Get (n=3) 1 2 2 1 2 20,1 3 0 1 2 2 0 O 1:4+4/4+4: l Fei et al. (2012) 
: f 57.2—69.7 62.7-71.5 d 
occidentalis (n=5) (n=5) = 3 2 3 2 2 - - - - 2 - 1 - - This study 
pingwuensis Ces a 0 2 3 1 2 3: L0,M 2 0 2 2 3 0 O 1:2+2/3+3: 1Fei et al. (2012) 


(n = 20) (n= 1) 
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species 
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(Continued Table 4) 


10 11 12 13 14 15 16 17 Sources 


46.9-51.3 50.8-53.9 


sikimmensis (n=8) (n=2) 1 1 - 1- - 3 LEM 
spinosus SE ER 3 SE 0 2 2 3 2 2 LI 
tuberculatus > are 3 "a 0 2 3 1 2 3 LI 

wanglangensis a pees l - - - -  - LU. UI) 
wuguanfui ey i 167 0 2 1.3 2 2 iti Ti 


disc ventrally set; lips with papillae, one uninterrupted 
row of keratodonts (denticles) and two interrupted rows of 
keratodonts on upper labium. Lower labium comprising 
one uninterrupted keratodont row followed by two 
interrupted keratodont rows. Tooth row formula (Dubois 
1995): 1:2+2/2+2:1 (Figure 7). Spiracle sinistral, single, 
midventral, posterodorsally directed; extended as a short 
tube, oval spiracular opening. Dextral vent tube opening 
at the margin of ventral fin. Intestinal coils visible. Tail 
musculature robust and greatly narrowing to tail tip; tail 
tip pointed; upper tail fin higher than lower tail fin, both 
fins converging at tip. 


Coloration in preservative 
The dorsal and lateral surfaces of the head, body and tail 
region are dark brown. The abdomen is dark brown and 
transparent, with clearly visible intestine coils; the ventral 
surface of tail creamy white. 


Comparisons 

Scutiger ghunsa sp. nov. is morphologically distinct from 
other Scutiger species found in Nepal and nearby regions 
by the following qualitative characters. The characters of 
Scutiger ghunsa sp. nov. are presented in parentheses; see 
also Table 4 and Appendix 1. 

Scutiger adungensis (Dubois, 1979), (distribution: 
only known from Adung valley of Myanmar) large sized 
(vs. medium sized), maxillary teeth present (vs. absent), 
toes weakly webbed (vs. absent), vocal sacs present (vs. 
absent). 

Scutiger bhutanensis (Delorme and Dubois, 2001), 
(distribution: Bhutan type locality unknown), large sized 
(vs. medium sized), maxillary teeth present (vs. absent), 
toes webbing weakly webbed (vs. absent), vocal sacs 


Jiang et al. (2012); 
Jiang et al. (2016); 
Subba et al. (2015) 


1 : 2+2 (3+3) / 
343:1 


1 
= 
N 
N 
U 
© 
© 


Jiang et al. (2012) 


1 : 2+2 (3+3) / This study; 
343): 1 Fei et al. (2012) 


W 
© 
N 
N 
w 
© 
© 


Fei et al. (2012) 


1 1:3+3/3+3: 1 Jiang et al. (2012) 


present (vs. absent), eye diameter larger than 20% of HL 
(vs. between 20—40% HL), interorbital distance shorter 
than IND (vs. longer), length of arm shorter than 20% of 
SVL (vs. between 20-40%), nuptial spines in breeding 
male on first two fingers (vs. presence in first three 
fingers), toes weakly webbed (vs. absent), vocal sacs 
present (vs. absent), pectoral and axillary glands similar 
in size (vs. two times larger). 

Scutiger boulengeri (Bedriaga, 1898), (Distribution: 
western Nepal, Sikkim, India, and southern Tibet to 
Sichuan, Yunnan, Gansu, and Qinghai, China, 3300-5270 
m elevation), female body size smaller (vs. medium size), 
snout shorter than 20% of HL (vs. between 20%—40% of 
HL), toes webbing weakly webbed (vs. absent), pectoral 
and axillary glands similar in size (vs. two times larger), 
spines on belly in breeding males (vs. absent), tadpole 


Figure 7 Mouth parts of tadpole of Scutiger ghunsa sp. nov. 
showing labial teeth. 
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Table 5 Measurements (mm) of the onomatophores of Scutiger ghunsa sp. nov. from Ghunsa, Taplejung, Nepal. 
Collection No. JRK- JRK- JRK- JRK- JRK- JRK- JRK- JRK- 
JRK2015-118 JRĶ2015-119 JRK2015-120 JRĶ2015-121 JRKĶK2015-122 JRĶ2015-123 JRK2015-124 JRK2015-125 
Holotype Paratype Paratype Paratype Paratype Paratype Paratype Paratype 
Male Male Male Male Male Female Female Female 
SVL 43.6 47.8 45.6 47.3 42.1 53.6 50.2 53.9 
HL 11.4 11.4 10.1 11.7 10.4 13.3 13.2 12.8 
HW 16.4 15.2 15.4 14.9 14.4 17.8 15.5 18.0 
SL 4.9 4.9 4.7 4.7 4.7 6.2 6.1 6.2 
ED 4.8 4.5 5.6 6.0 4.8 7.2 6.2 5.0 
ENL 3.1 37 2.3 2.4 2.7 3.2 3.4 3.5 
UEW 44 2.4 2T 4.4 4.4 4.6 4.3 4.6 
IOW 5.4 4.7 4.3 5.3 6.5 5.0 5.6 4.0 
IND 4.1 3.7 3.5 3.5 3.6 4.0 4.9 4.6 
FAL 10.0 11.5 9.4 12.1 8.9 10.1 11.9 12.0 
HAL 23.4 21.8 21.4 19.9 20.3 25.1 21.6 23.7 
Fl 3.1 4.6 4.7 3.9 4.0 4.1 4.7 4.4 
F2 3.6 3.7 6.1 5.0 4.4 4.5 4.7 5.4 
F3 5.9 6.1 7.4 6.9 TS 8.3 8.3 10.8 
F4 4.8 5.4 4.5 6.1 5.8 7.8 6.1 7.6 
THL 17.5 19.4 19.4 16.5 13.4 20.6 16.3 20.0 
SHL 17.1 17.0 17.6 13.6 16.0 17.4 15.3 18.7 
TL 11.1 1321 11.4 8.1 10.4 13.1 12.3 7.8 
FOL 19.3 21.4 20.1 18.9 18.6 22.4 22.2 22.8 
Tl 3.7 4.1 7.0 3.6 4.0 4.5 6.6 3.6 
T2 5.4 7.8 8.1 4.1 5.2 6.7 7.7 5.4 
T3 8.6 10.3 TL 5.8 8.3 9.6 10.2 9.3 
T4 11.5 15.9 15.6 12.0 9.2 12.4 13.6 12.9 
T5 6.9 10.4 E 5.6 8.9 9.1 77 9.7 


keratodont formula (1:5+5(6+6)/5+5(6+6):1). 

Scutiger brevipes (Liu, 1950), (distribution: Taining, 
Daofu County, Sichuan, China), large sized (vs. medium 
sized), nuptial spines on breeding male on first two fingers 
(vs. present on first three fingers), toes weakly webbed 
(vs. absent). 

Scutiger chintingensis (Liu and Hu, 1960), 
(distribution: Emei Mountain region, Hongya and 
Wenchuan counties, Sichuan, China), maxillary teeth 
present (vs. absent), arm shorter than 20% of SVL (vs. 
between 20-40%), toes webbed (vs. webbing absent), 
tadpole keratodont formula (1:3+3/2+2:1). 

Scutiger glandulatus (Liu, 1950) (distribution: 
Northern Yunnan, western Sichuan and southern Gansu, 
China.), large sized (vs. medium sized), nuptial spines in 
breeding male on first two fingers (vs. presence on first 
three fingers), presence subarticular tubercles (vs. absent), 
toes weakly webbed (vs. absent), tadpole keratodont 


formula (1:4+4/4+4:1). 

Scutiger gongshanensis (Yang and Su, 1979) 
(distribution: Gongshan and Biluoxshan in Yunnan, 
China), presence of maxillary teeth (vs. absent), nuptial 
spines in breeding males on first two fingers (vs. presence 
on first three fingers), a single pectoral gland in breeding 
male (vs. two), tadpole keratodont formula (1:3+3/3+3:1). 

Scutiger jiulongensis (Fei, Ye, Cand Jiang, 1995) 
(distribution: Tanggu, Juilong County, Sichuan, China), 
large sized (vs. medium sized), nuptial spines in breeding 
males on first two fingers (vs. presence on first three 
fingers), subarticular tubercles present (vs. absent), toes 
webbed (vs. absent), tadpole keratodont formula (1:3+3/ 
4+4(3+3):1). 

Scutiger liupanensis (Huang, 1985), (distribution: 
Liupanshan, Jingyuan County, Ningxiz Huiza Zizhiqu 
Autonomous Region, China), presence of maxillary teeth 
(vs. absent), toes weakly webbed (vs. absent), pectoral 
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and axillary glands similar in size (vs. two times larger), 
spines on belly in breeding males (vs. absent), tadpole 
keratodont formula (1:6+6/6+6:1). 

Scutiger maculatus (Liu, 1950), (distribution: Garza 
County, Sichuan, and Jiangda, Tibet, China), large body 
size (vs. medium size), maxillary teeth present (vs. 
absent), toes weakly webbed (vs. absent), pectoral and 
axillary glands similar in size (vs. two times larger). 

Scutiger mammatus (Günther, 1896), (distribution: 
Eastern Tibet, Qinghai and Sichuan, China), body large 
sized (vs. medium sized), nuptial spines in breeding males 
on first two fingers (vs. presence on first three fingers), 
presence of subarticular tubercles (vs. absent), toes 
webbed (vs. absent), a single pectoral gland in breeding 
male (vs. two), tadpole keratodont formula (1:4+4/5+5 
(4+4):1). 

Scutiger muliensis (Fei and Ye, 1986), (distribution: 
Muli County, Sichuan, China), body large sized (vs. 
medium sized), nuptial spines in breeding males on 
first two fingers (vs. present on first three fingers), toes 
weakly webbed (vs. absent), a single pectoral gland in 
breeding males (vs. two), tadpole keratodont formula 
(1:4+4/4+4:1). 

Scutiger nepalensis (Dubois, 1974), (distribution: 
Khaptad, western Nepal), body large sized (vs. medium 
sized), eye diameter larger than 20% of HL (vs. between 
20—40% HL), interorbital distance shorter than IND (vs. 
longer), length of arm shorter than 20% of SVL (vs. larger 
than 40% of SVL), pectoral and axillary glands similar in 
size (vs. pectoral gland two times larger). 

Scutiger ningshanensis (Fang, 1985), (distribution: 
Ningshan County in southern Shaanxi Province, China), 
maxillary teeth present (vs. absent), toes weakly webbed 
(vs. absent), pectoral and axillary glands similar in size 
(vs. pectoral gland two time larger), spines on belly in 
breeding males (vs. absent), tadpole keratodont formula 
(1:5+5/5+5:1) 

Scutiger nyingchiensis (Fei, 1977), (distribution: 
Western Nepal, Kashmir, Ladakh, and Tibet, China), 
large body size (vs. medium size), presence of maxillary 
teeth (vs. absent), eye diameter larger than 20% of HL 
(vs. between 20—40% of HL), toes webbed (vs. absent), 
tadpole keratodont formula (1:4+4/4+4:1). 

Scutiger occidentalis (Dubois, 1978) (distribution: 
Jammu Kashmir, India), large body size (vs. medium 
size), snout shorter than 20% of HL (vs. between 20—40% 
of HL), eye diameter larger than 20% of HL (vs. between 
20—40% of HL), interorbital distance shorter than IND 
(vs. longer), length of hand shorter than 20% of SVL (vs. 
between 20—40% of SVL). 
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Scutiger pingwuensis (Liu, Hu, and Tian, 1978), 
(distribution: Pingwu County, Sichuan, China), large body 
size (vs. medium size), length of inner metatarsal tubercle 
equal to first toe length (vs. less than 50% or shorter), toes 
weakly webbed (vs. absent), tadpole keratodont formula 
(1:5+5/3+3:1). 

Scutiger sikimmensis (Blyth, 1854), (distribution: 
Sikkim and Meghalaya (India), Bhutan and southern 
Tibet, China), maxillary teeth present (vs. absent), snout 
length shorter than 20% of HL (vs. between 20—40% ), 
subarticular tubercles present (vs. absent), toes weakly 
webbed (vs. absent), tadpole keratodont formula 
(1:2+2(3+3)/3+3:1). 

Scutiger spinosus (Jiang, Wang, Zou, Yan, Li and Che, 
2016) (distribution: Medog, Tibet, China), body large 
sized (vs. medium sized), eye diameter larger than 20% 
of HL (vs. between —40% of HL), interorbital distance 
shorter than IND (vs. longer), length of hand between 
20-40% of SVL (vs. longer than 40% of SVL), inner 
metatarsal tubercle equal length of first toe (vs. less than 
50% or shorter), toes weakly webbed (vs. absent). 

Scutiger tuberculatus (Liu et al., 1979), (distribution: 
Southern Sichuan and north-central Yunnan, China), 
body large sized (vs. medium sized), nuptial spines in 
breeding males on first two fingers (vs. presence on first 
three fingers), toes weakly webbed (vs. absent), tadpole 
keratodont formula (1: 2+2(3+3)/3+3:1). 

Scutiger wanglangensis (Ye and Fei, 2007), 
(distribution: Pingwu and Nanping counties, Sichuan, 
China), body large sized (vs. medium sized), toes weakly 
webbed (vs. absent), pectoral and axillary glands similar 
in size (vs. pectoral gland two times larger), spines on 
belly in breeding males (vs. absent). 

Scutiger wuguanfui (Jiang, Rao, Yuan, Li, Hou, Chen 
and Che, 2012) (distribution: Medog, Muotuo County, 
Xizang, China), body large sized (vs. medium sized), eye 
diameter larger than 20% of HL (vs. between 20%—40% 
of HL), interorbital distance shorter than IND (vs. longer), 
length of hand between 20%—40% of SVL (vs. longer 
than 40% of SVL), inner metatarsal tubercle equal to 
length of first toe (vs. less than 50% of first toe or shorter), 
toes weakly webbed (vs. absent), vocal sacs present (vs. 
absent), tadpole keratodont formula (1:3+3/3+3:1). 


Distribution and natural history 

Scutiger ghunsa sp. nov. is currently known only from 
the type locality (Figure 8) and could not be found at 
several localities from Makalu Barun National Park 
and Kanchenjunga Conservation Area in the eastern 
Nepal. The new species was recorded from a mixed 
rhododendron forest with low canopy cover with a slow- 
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running stream and stagnant water pools. We encountered 
only five specimens at night, hiding on boulders and in 
rock crevices. We found an adult male with clutches on its 
back hidden in the rock crevices. Some frogs are known 
to guard egg clutches (Zheng et al., 2011), however, 
parental care in the genus Scutiger has not been reported 
so far. But, we have no further evidence to confirm if 
this behaviour was accidental or if it can be considered 
parental care. 

Tadpoles were free moving and were highly 
photosensitive at night. They usually take refuge in rock 
crevices and vegetation when a torch light is pointed at 
them. The water was clear, and soil was acidic (pH 5.4) 
in the habitat. No calls were heard. Type locality lies 
inside the Kanchenjunga Conservation Area, in a semi- 
forested part, near a human trail and Ghunsa Khola. 
Therefore, considering the lack of comprehensive records 
of populations of the new species, we suggest that the 
Red List category of Scutiger ghunsa sp. nov. is Data 
Deficient according to IUCN (2001) guidelines. 


4. Discussion 


More than 57% (12 out of 21) of known species of 
Scutiger occur at the edge of the Sichuan basin in western 
China and many of them are reported only from their 
type-locality (Fei et al., 2012). Recently, the genus 
Scutiger has been divided in to three major groups: the 
greater Himalayan, Tsinling Mountaing and Sichuan 
Basin, and the Tibet and Hengduan Mountain regions 
based on molecular phylogenetic analysis (Hofmann 
et al., 2017). Scutiger species that are distributed on the 
southern slope of the Himalaya belongs to the greater 
Himalayan species (S. nepalensis, S. sikkimensis, S. 
occidentalis, and S. ghunsa sp. nov.). Our molecular 
results suggested that S. nepalensis is the closest congener 
to the new species; this species is distributed at elevations 
of 2700-4500 in western Nepal. Scutiger ghunsa sp. nov 
showed substantial differentiation from S. nepalensis in 
both morphological and molecular analysis. Among the 
known Scutiger species, S. sikkimensis occurs in closest 
proximity to the new species and is distributed in the 
Dudh Koshi and Arun rivers approximately 40-60 km 
west from the type locality of S. ghunsa sp. nov. However, 
their ranges are greatly separated from each other by 
steep river gorges and mountains. This study identified 
three distinct divergent lineages of S. sikkimensis, 
two of them are from Dhudh Koshi and Arun River, 
Nepal and one from Yadong, Tibet, China. There was 
substantially low genetic divergence in 16S gene between 


Figure 8 A. Habitat and B. microhabitat of Scutiger ghunsa sp. 
nov. at Kanchenjunga Conservation Area, Taplejung district, Ghunsa 
valley, Nepal. 


S. sikkimensis from Sikkim, India and Dhudh Koshi and 
Arun River, Nepal. Due to high genetic divergence that 
exists in this species group, further study is required on 
these two populations from Yadong, China and Surkie La, 
Dudhkoshi River (west) to solve the current taxonomic 
ambiguities. 

This study found that S. boulengeri is distributed 
widely to the alpine desert across the Tibetan plateau 
of China, India (Sikkim) and Nepal (Mustang). 
However, previous studies suggested that S. boulengeri 
is a paraphyletic species complex with a high rate of 
mitochondrial gene introgression following hybridization 
(Chen et al., 2009). The molecular phylogeny of the genus 
Scutiger indicated that this group appears to have highly 
conserved morphological characters and suggests the 
presence of many many cryptic species (Hofmann er al., 
2017). Of the currently recognized species of Scutiger, S. 
bhutanensis and S. adungensis are geographically closer 
to Scutiger ghunsa sp. nov. Scutiger bhutanensis is known 
from the imprecise type locality from Bhutan, without any 
habitat and other geographical characteristics (Delorme 
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and Dubois 2001). The second species, S$. adungensis, 
from Adung Valley of Myanmar near to China, boarder at 
3,650 m, is also reported from a single locality. 

Previously, four species of Scutiger were reported 
from Nepal based on external morphology (Dubois 1974; 
1987; Nanhoe and Ouboter 1987). Past studies were 
limited by several factors such as limited sampling in 
restricted geographic areas or limited species included in 
the morphological comparisons. In recent years, several 
cryptic taxa have been identified based on molecular 
methods. For example, Hofmann er al., (2017) reported 
highly divergent lineages of S. nepalensis and S. 
sikkimensis from western and eastern Nepal Himalayas 
that might be a potential new species. In the framework 
of an integrative taxonomic approach (Dayrat 2005; 
Schlick-Steiner et al., 2010) to describe biodiversity 
and understand the diversification of Scutiger across 
the Himalayan range, not only is further fieldwork 
needed throughout the range, but diverse approaches 
(morphological studies including historical specimens, 
in particular types, behavioral studies, specially analyses 
of mating calls and mating behavior, or cytology) should 
be used in combination with molecular studies. In recent 
publications of the amphibians of the eastern Himalayas 
of Nepal, two new species of amphibians, Tylototriton 
himalayanus (Khatiwada et al., 2015) and Microhyla 
taraiensis (Khatiwada et al., 2017) were identified, 
indicating the presence of several unrecognized cryptic 
species. With the addition of Scutiger ghunsa sp. nov. 
the genus Scutiger now contains 24 species with 5 in 
Nepal. The description of a new species highlights the 
still poor knowledge of amphibians in the area and points 
to the possibility of the existence of more cryptic taxa. 
Because Scutiger ghunsa sp. nov. is currently known only 
from type locality, efforts are needed to find additional 
populations of the new species and assess its conservation 
status. 
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Appendix 


Table S1 Currently recognized species of genus Scutiger and known distributional ranges. 


Species Author and date Type locality Known distribution Reference 
Scutiger : Adung valley of Myanmar, near : : 
dancers Dubois, 1979 Myanmar and China border Myanmar and China border Dubois (1979) 
Scutiger Delorme and Dubois, Delorme and Dubois 
bhutanensis 2001 pintan Bhutan (2001) 

P e North-western Nepal, Sikkim (India), and Qinghai Fei et al. (2012); 
ae Bedriaga, 1898 a EEN Tibetan plateau to Sichuan, and Gansu Provinces of Frost (2018), 

g anaes China, 3300-5270m. This study 
Scutiger ; Taining, Daofu County, e ‘ e Fei et al. (2012); 
brevipes eee Sichuan, China Sichuan Bovine, Cima Frost (2018) 
Scutiger i Chinting, Mt. Omei, Sichuan, : S z Fei et al. (2012); 
chintingensis Piven Hu, 1300 China, 3050m altitude SE Frost (2018) 
Scutiger > Hopachai, Nibashan, Sichuan, Normer Yunnan (zhongdian) through westërn and Fei et al. (2012); 

Jlandulatis Liu, 1950 China, 2590m altitude central Sichuan to extreme southern Gansu (Wenxian), p st (2018) 
S i China, 2200 to 4000m. S 
Scutiger Yane and Su. 1979 12th Bridge, Gongshan Xian, Gongshan and Biluoxueshan in northwestern Yunnan, Fei er al. (2012); 
gongshanensis 8 e Yunnan, China, 2750m altitude China and Myanmar, 2500-3850 m. Frost (2018) 
Scutiger : ; Juilong, Kangding, Sichuan, d e i Fei et al. (2012); 
jiulongensis Fei, Ye, and Jiang, 1995 China Tanggu, Juilong County, Sichuan, China, 3210m. Frost (2018) 
Scutiger Mt Tiupan, Ningxia Hur Liupan, Jingyuan, Ningxia, Pingliang and Zhuanglang Fei et al. (2012); 
i s Huang, 1985 Autonomous Region), China $ : 
liupanensis . counties, Gansu, China, 1900—2500m. Frost (2018) 
2100m altitude 
Scutiger Liu. 1950 Hasa, Kantze, Sichuan, China, Garze County in northwestern Sichuan, and Jomda, in Fei et al. (2012); 
maculatus > 3350m altitude eastern Tibet, China, 3300—3500m. Frost (2018) 
Eastern Tibet (Baxoi, Qamdo, Jomda, Markam, 
i Riwoge, and Zayu counties), Qinghai (Baima, ; i 
geneeë é Giinther, 1896 Tungsolo, Sichuan, China Nanggen, and Yushu counties), Sichuan (Dege, a a Se ae 
j Batang, Kangding, and Luhuo counties), and Yunnan 3 
(Deqen County), China. 
Scutiger Fei and Ye. 1986 Muli County, Sichuan, 3200m Muli County, southwestern Sichuan, China, 3050- Fei et al. (2012); 
muliensis $ g altitude 3400m. Frost (2018) 
Scutiger : Khaptad, Chainpur, Nepal o Dubois (1974); 
Dë? Dubois, 1974 2950m altitude Kanjiroba, Saipal Mountain, Jumla western Nepal Hofmann et al. (2017) 
Scutiger Fane, 1985 Xunyangba, Ningshan County, Ningshan County in southern Shaanxi, east to southernFei et al. (2012); 
ningshanensis & Shaanxi, China, 2,550m altitudeHenan, China, 1970—2550m. Frost (2018) 
Scutiger , Norecht Tibet China; 3040m O eeter Nepal anid southeastern Tibet (Bomi. ` et aJ (ION: 
ee E Fei, 1977 altitude Lhozhag, Mainling, Medog, Nangxian, Nyingchi, and Frost (2018) 
ý Yadong counties), China, 2730-4560m. 

3 Shukdhari, (Jammu and Gg : ; : Hofmann et al. 
Scutiger Dubois, 1978 Kashmir, India, 2920-2940m Gilgit Baltistan, Pakistan, and Ladakh in Kashmir, (2017); 
occidentalis j India. j 

altitude This study 
Scutiger Liu and Tian, 1978 Wang-ba-chu, Pingwu County, Wenxian County, southern Gansu, and Pingwu and Fei et al. (2012); 
pingwuensis S > Sichuan, China, 2200m altitude Beichuan counties, Sichuan, China, 2100—2200m Frost (2018) 

! PEP : Fei et al. (2012); 
Scutiger aah : Eastern Nepal, Sikkim and Meghalaya, India, Bhutan ae 
sikimmensis Fe Sikkim, India and adjacent southern Tibet, China, 2860-4116m Frost (2018); 

Jiang et al. (2016) 
Scutiger Jiang, Wang, Li, and 62K, Medog, Tibet, China, High altitude forest around Medog, Sela Pass, adjacent 


spinosus 


Che, 2016 


2705m altitude 


southern Tibet, China, 2705-4169m 


Jiang et al. (2016) 


(Continued Table S1) 


Species Author and date Type locality Known distribution Reference 
Scutiger Liu and Fei, 1979 Pu-xiong of Yuexi county, Mianning, Xichang, Yanbian, Yuexi, and Zhaojue Fei et al. (2012); 
tuberculatus í Sichuan, altitude 3000m, China.counties, Sichuan, 2600-3750m Frost (2018) 


Scutiger Ye and Fei, 2007 


Wenxian County in Gansu Province, and Pingwu and 


Wangbachu of Pingwu County, Nanping counties, Sichuan Province, China, 2200- 


Fei et al. (2012); 


wanglangensis Sichuan, China, 2500m altitude 2800m Frost (2018) 
Scutiger Jiang, Rao, Yuan, Wang,62K of Medog, Medog County, 
Wee fui Li, Hou, Che, and Che, Tibet Autonomus Region, Medog County in Tibet of China Jiang et al. (2012) 


2012 


China, 2705m altitude 


Table S2 Specimens used for morphological comparisons. JRK refers to collection number of Janak Raj Khatiwda, CIB refers to Chengdu 
Institute of Biology, MMHN to Muséum national d’histoire naturelle, Paris, France and NMBS to Natural History Museum of Basel, 
Switzerland. 


Museum collection number 


Species name Locality 
Male Female 
JRK2015-193 
JRK2015-194 JRK2015-198 
S. ghunsa sp. nov. Ghunsa, Taplejung, Nepal JRK2015-195 JRK2015-199 
JRK2015-196 JRK2015-200 
JRK2015-197 BRG 
J 16-11 
S. boulgeneri Muktinath, Mustang, Nepal E JRK2016-12 
JRK2016-13 
CIB790297 
CIB790286 CIB790296 
Qinghai, China CIB790287 CIB790573 
CIB790288 CIB790297 
CIB790285 
CIB25933 
CIB25939 CIB25936 
S. nyingchiensis Bomi, Tibet, China CIB25937 CIB125942 
CIB25934 CIB25947 
CIB15946 
CIB83434 
CIB83435 
S. pingwuensis Pingwu, Sichuan, China CIB83431 
CIB84467 
CIB84468 
ee CIB87182 C1887183 
S. tuberculatus Yaan, Sichuan China CIB87181 
CIB26153 CIB87179 
CIB26154 
CIB24672 
CIB24676 CIB83442 
S. chintingensis Wawushan Sichuan, China CIB24675 CIB95149 
CIB24674 CIB95146 
CIB24673 
CIB84861 
CIB84860 CIB26051 
S. mammatus Batang, Sichuan, China CIB84862 CIB26048 
CIB84863 CIB26036 
CIB84853 
CIB94965 
CIB87177 CIB94969 
S. glandulatus Yaan, Sichuan China CIB87173 CIB94967 
CIB25867 CIB25870 
CIB25866 
MNHN-RA-1974.1096 
Š nepalensis Khaptad, Nepal MNHN-RA-1989.3361 MNHN-RA-1974.1098° 
` MNHN-RA-1989.3362 MNHN-RA-1974.1097 
MNHN-RA-1974.1095 
MNHN-RA-1977.1097 MNHN-RA-1977.1120 
MNHN-RA-1977.1080 MNHN-RA-1977.1114 
S. occidentalis Jammu and Kashmir, India AE EEN 
i MNHN-RA-1977.1105 MNHN-RA-1977.1121 
MNHN-RA-1977.1089 MNHN-RA-1977.1116 
MNHN-RA-1977.1069 MNHN-RA-1977.1118 
S. bhutanesis Bhutan NMBS09147 


Note: Juveniles 


va Scutiger ghunsa sp. nov. (JRK193) Ghunsa, Taplejung, Nepal 
Scutiger ghunsa sp. nov. (JRK194) Ghunsa, Taplejung, Nepal 
Scutiger ghunsa sp. nov. (TRK195) Ghunsa, Taplejung, Nepal 
— Scutiger ghunsa sp. nov. (JRK197) Ghunsa, Taplejung, Nepal 
- Scutiger ghunsa sp. nov. (JRK199) Ghunsa, Taplejung. Nepal 
Scutiger ghunsa sp. nov. (JRK206) Ghunsa, Taplejung, Nepal 
Scutiger ghunsa sp. nov. (JRK202) Ghunsa. Taplejung, Nepal 
` Scutiger ghunsa sp. nov. (JRK203) Ghunsa. Taplejung. Nepal 
Scutiger ghunsa sp. nov. (JRK196) Ghunsa, Taplejung, Nepal 
- Scutiger ghunsa sp. nov. (JRK198) Ghunsa, Taplejung, Nepal 
Scutiger ghunsa sp. nov. (JRK200) Ghunsa, Taplejung, Nepal 
Scutiger ghunsa sp. nov. (JRK201) Ghunsa, Taplejung, Nepal 
Scutiger ghunsa sp. nov. (JRK204) Ghunsa, Taplejung, Nepal 
Scutiger ghunsa sp. nov. (JRK205) Ghunsa, Taplejung, Nepal 
~ Scutiger ghunsa sp. nov. (JRK207) Ghunsa, Taplejung, Nepal 


0.68 0.99 . nepalensis, Kanjiroba, Jumla, Nepal 

. nepalensis, Saipal-Juphal, Nepal 

. nepalensis, Dhauli lake, Jumla, Nepal 

. nepalensis, Dhaulung-North, Jumla, Nepal 
. nepalensis, Dhaulung-South, Jumla, Nepal 


S 
S 
S 
S 
S 
—— S. nepalensis, Saipal-chainpur, Jumla, Nepal 
S 
S 
S 
S 
S 


, nyingchiensis, Nyingchi, Tibet, China 
` S. gongshanensis, Gongshan, Yunnan, China 
. spinosus, Medog, Nyingchi, Tibet, China 
S. boulengeri (JRK211), Muktinath, Mustang, Nepal 
. boulengeri (JRK212), Muktinath, Mustang, Nepal 


S. boulengeri (JRK214), Muktinath, Mustang, Nepal 
S. boulengeri (JRK215), Muktinath, Mustang, Nepal 
. boulengeri, Tageija, Tibet, China 


S 
S 
0.60 0.64 S 
S. boulengeri (JRK213), Muktinath, Mustang, Nepal 
S 
S, 
S 
S 


S. boulengeri, Nanqian, Zhaduo, Quinghai, China 
S. boulengeri, Suojia, Zhiduo, Qinghai, China 

. boulengeri, Suojia, Zhiduo, Qinghai, China 

S. boulengeri, Damxung, Tibet, China 

. boulengeri, Lasha, Tibet, China 


0:63 S. glandulatus, Garze, Sichuan, China 


S. jiulongensis, Ganzi, Sichuan, China 


S 
S 
S. 
S. liupanensis, Liupan Mountain, Gansu, China 
S. occidentalis, Deosai Plains, Pakistan 
S. sikkimensis. Nepal, Arun, Kongma Danda, Nepal 
~~ S. sikkimensis. Taksindo, Dudhkoshi river-East, Nepal 
S. sikkimensis, Surkie La, Dudhkoshi river-West, Nepal 
1 L~~ S, chintingensis, Tianquan, Sichuan, China 
~ S. sikkimensis, Yadong, Tibet, China 
, ningshanensis, Taibai, Shaanxi, China 
. wuguanfui, Medog, Nyingchi, Tibet, China 
1-——_----------- -—-——- Megophrys jingdongensis, Lao Cai, Sapa, Vietnam 
1 TU -—-—- Megophrys omeimontis, Emei Shan, Sichuan, China 
Megophrys maosonensis, Xiaogiaogou Nature Reserve, Yunnan, China 
— Duttaphrynus melanostictus, India 


0.07 


Figure S1 Bayesian inference (BI) tree based on COI sequences. Values on branches of the tree are Bayesian posterior probabilities (bpp) 
followed by ML bootstrap support (bs) below. Asterisk (^) denotes bpp > 90. 


S. boulengeri, Zhaduo, Qinghai, China 
-- S. boulengeri, Zhiduo, Qinghai, China 
- S. boulengeri, Sikkim, India 
- S. boulengeri, Damxung, Tibet, China 
- S. boulengeri, Muktinath, Mustang, Nepal 
- S. boulengeri, Lasha, Tibet, China 
` S. boulengeri, Zhiduo, Qinghai, China 
` S. boulengeri, Tageija, Tibet, China 
- S. boulengeri, Muktinath, Mustang, Nepal 
- S. boulengeri, Muktinath, Mustang, Nepal 
S. boulengeri, Muktinath, Mustang, Nepal 
~ S. boulengeri, Muktinath, Mustang, Nepal 
- S. liupanensis, Liupan Mountain, Gansu, China 
. muliensis, Yanyuan, Sichuan, China 
. tuberculatus, Xichang, Sichuan, China 
. mammatus, Baiyu, Sichuan, China 
~ S. glandulatus, Garze, Sichuan, China 
ES S. wuguanfui, Medog, Nyingchi, Tibet, China 

30 S. ningshanensis, Taibai, Shaanxi, China 

EE S. chintingensis, Tianquan, Sichuan, China 
- S. occidentalis, Deosai Plains, Pakistan 
` Scutiger ghunsa sp. nov. (JRK196) Ghunsa, Taplejung, Nepal 
` Scutiger ghunsa sp. nov. (IRK198) Ghunsa, Taplejung, Nepal 
` Scutiger ghunsa sp. nov. (JRK200) Ghunsa, Taplejung, Nepal 
~ Scutiger ghunsa sp. nov. (IRK201) Ghunsa, Taplejung, Nepal 
~ Scutiger ghunsa sp. nov. (JRK202) Ghunsa, Taplejung, Nepal 
` Scutiger ghunsa sp. nov. (JRK203) Ghunsa, Taplejung, Nepal 
` Scutiger ghunsa sp. nov. (IRK204) Ghunsa, Taplejung, Nepal 
` Scutiger ghunsa sp. nov. (JRK205) Ghunsa, Taplejung, Nepal 
` Scutiger ghunsa sp. nov. (JRK207) Ghunsa, Taplejung, Nepal 
` Scutiger ghunsa sp. nov. (JRK193) Ghunsa, Taplejung, Nepal 
` Scutiger ghunsa sp. nov. (IRK194) Ghunsa, Taplejung, Nepal 
` Scutiger ghunsa sp. nov. (IRK195) Ghunsa, Taplejung, Nepal 
` Scutiger ghunsa sp. nov. (JIRK197) Ghunsa, Taplejung, Nepal 
` Scutiger ghunsa sp. nov. (JRK199) Ghunsa, Taplejung, Nepal 
~ Scutiger ghunsa sp. nov. (JRK206) Ghunsa, Taplejung, Nepal 
5. nepalensis, Saipal-Chainpur, Jumla, Nepal 
. nepalensis, Saipal-Juphal, Nepal 
. nepalensis, Kanjiroba, Jumla, Nepal 
. nepalensis, Dhauli lake, Jumla, Nepal 
. nepalensis, Dhaulung-North, Dolpa, Nepal 
. nepalensis, Dhaulung-South, Dolpa, Nepal 
. sikkimensis, Surkie La, Dudhkoshi river-West, Nepal 
. sikkimensis, Arun, Kongma Danda, Nepal 
. sikkimensis, Sikkim, India 
. sikkimensis, Taksindo, Dudhkoshi river-east, Nepal 
-~ O. jingdongensis, Sichuan, China 
" O. major, Sichuan, China 
-~ O. xiangchengensis, Yunnan, China 
~~ O. rhodostigmatus, Guizhou, China 
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57 
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81 
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37 


91 


Figure S2 Maximun likelihood tree based on 16S gene sequences. Values on branches of the tree are ML bootstrap supports (bs). 


